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Abstract 

We have investigated the electronic structure of the zigzag ladder (chain) compound SrCu02 
combining polarized optical absorption, reflection, photoreflectance and pseudo-dielectric function 
measurements with the model calculations. These measurements yield an energy gap of 1.42 eV 
(1.77 eV) at 300 K along (perpendicular) to the Cu-0 chains. We have found that the lowest 
energy gap, the correlation gap, is temperature independent. The electronic structure of this oxide 
is calculated using both the local-spin-density-approximation with gradient correction method, and 
the tight-binding theory for the correlated electrons. The calculated density of electronic states 
for non-correlated and correlated electrons shows quasi-one-dimensional character. The correlation 
gap values of 1.42 eV (indirect transition) and 1.88 eV (direct transition) have been calculated with 
the electron hopping parameters t = 0.30 eV (along a chain), t yz = 0.12 eV (between chains) and 
the Anderson-Hubbard repulsion on copper sites U= 2.0 eV. We concluded that SrCu02 belongs 
to the correlated-gap insulators. 

PACS numbers: 78.40.-q; 71.27.+a; 71.15.Ap; 71.20.Ps; 
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I. INTRODUCTION 



Strontium copper oxide, SrCu02, belongs to the new family of quasi-one- dimensional 
(ID) insulators whose properties have been a subject of intensive studies in recent years 
[| f|, ||, H |7|, |8|, [| [10], |TT|, [T2|]. This oxide, grown in single crystalline form under ambient 
pressure, has an orthorhombic unit cell (space group D^ h - Cmcm) with parameters a=3.577 
A, b=16.342 A, c= 3.9182 A, Z=4 [TJ. SrCu0 2 has a unique structure consisting of 
CuC>4 squares, mutually connected via common edges, that form double copper zigzag chains 
(Fig. 1). 

Magnetic susceptibility measurements of SrCuC>2 have revealed that the Cu 2+ moments 
order antiferromagnetically below ~ 2K ||. The exchange interaction energy J is estimated 
to be 2100 ± 200 K ||. Very recently, a static disorder spin structure (spin freezing) rather 
than static three dimensional long range spin order is found in SrCuC>2 using neutron spec- 
troscopy An angle-resolved photoemission spectroscopy (ARPES) study of SrCuC>2 
shows the spin-charge separation in this oxide, as a consequence of electron correlations 
||. [|. The optical phonons in SrCu0 2 have been investigated by measuring Raman scatter- 
ing H |7|, ||, H and far-infrared reflectivity spectra on policrystalline |7j and single crystal || 
samples. 

Electronic structure of SrCuC>2 has been calculated using the linear-augmented-plane- 
wave (LAPW) method within the local-density approximation (LDA) to the density func- 
tional theory (10| . As noted in Ref . [|10j for this type of calculations the agreement with the 
ARPES data takes place only in the region far from the Fermi level. Most recent band 
structure calculations [pi] , |I~2f by the local-spin-density-approximation (LSDA), including 
an on-site Coulomb repulsion (LSDA+U method) have produced an insulating gap of 1.63 
eV for U=12.3 eV |ll| and 2.36 eV for U=5 eV fl2| . Due to the lack of experimental data 
for the energy gap of SrCu02, it was not possible to compare these data with experiments. 



Quite recently, Rosner et al. |1| calculated the electronic structure of SrCu02 using the 
LSDA method, also. They found that the gap in the antiferromagnetic phase of SrCu02, 
reported in previous publications [|TT], [12| , is probably artificial and appears due to too small 
number of the k-points in the Brillouine zone taken in the computation procedure. 

In this paper we have applied complementary optical spectroscopy techniques to study 
the electronic structure of SrCu0 2 . We have obtained an energy gap of 1.42 eV (1.77 eV) 
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at 300 K along (perpendicular to) the Cu-0 chains. These values are compared with the 
results of the electronic structure calculations completed by us. We have established that 
the tight-binding method for the correlated electrons, with the hopping parameters t = 
0.30 eV and t yz = 0.12 eV along legs and between legs, respectively, (see Fig. 1) and the 
Anderson-Hubbard parameter U = 2.0 eV yields the energy gap values in an agreement 
with the experimental results. We also have determined the density of electronic states for 
non-correlated and correlated electrons, which show the one-dimensional nature. According 
to the electronic structure calculations and the experimental findings, we concluded that 
SrCuC>2 is a correlated-gap insulator. 

The remainder of this paper is organized as follows. Experimental details are given in 
Sec. II. The results of electronic structure calculations for non-correlated electrons (using the 
LSDA with gradient corrections to density functional theory) are described in Sec. III. A. The 
reason why we use this method is to obtain the data necessary for the determination of the 
correlated electronic structure and to compare our results with the previously published ones 



T0| , O , [L2L [16|]. In Sec. III.B, we applied the tight-binding theory for correlated electrons 



taking into account the realistic crystalline structure of SrCu02- We got analytically the 
energy dispersion relations, the density of states of correlated and non correlated electrons 
as well as the correlation gap. Experimental measurements of dielectric function, reflectivity, 
photoreflectance and optical transmission spectra are given in Sec. IV. Section V contains 
the discussion of experimental and theoretical results and the conclusions. 



II. EXPERIMENTAL DETAILS 



The present work was performed on a single crystal sample with a size of about 15, 2 and 
6 mm along the a, b and c axes, respectively. We used several optical spectroscopy tech- 
niques. The pseudo-dielectric function was measured with the help of a rotating-polarizer 
(analyzer) ellipsometer. A Xe-lamp was used as a light source, a double monochromator 
with 1200 lines/mm gratings and an S20 photomultiplier tube as a detector. The polarizer 
and analyzer were Rochon prisms. The measurements were carried out in the 1.2-5.6 eV 
energy range. For the energies below 1.6 eV we used a halogen lamp as a light source and 
Si-photo diode as a detector. Optical reflectivity and transmission spectra were measured at 
room and liquid helium temperature in the 200-2500 nm spectral range with Perkin - Elmer 
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Lambda 19 spectrophotometer. Photoreflectance measurements were performed at room 
temperature under incident polarized light parallel to the a- and the c-axes. The reflected 
light was dispersed through a single 1/2 m monochromator and sinchronously detected by 
Si-photodiode. 



III. ELECTRONIC STRUCTURE CALCULATIONS 



A. LSDA with gradient corrections 

In the CuCVlayered oxides, the electron dispersion branches at the Fermi level are de- 
rived from the 3d x 2_ y 2 copper orbitals, mixed with the corresponding O 2p orbitals. Other 
occupied 3d Cu orbitals (d 3z 2_ r 2,d yz ,d xz ,d xy ) are located below the Fermi level and they 
are strongly hybridized with oxygen bands. These cuprates are known as charge-transfer 
insulating oxides. Contrary to the CuCvlayer cuprates (tetragonal crystal structure ox- 
ides), SrCu02 has an orthorhombic crystal structure with the Cu-0 zigzag chains and, 
consequently, with the different electronic structure, as it will be shown later. 

For the ab-initio calculation of the electronic structure of SrCu02 we used WIEN97 
software package |17]]. The program allows us to compute the electronic structure within 



the density functional theory by applying the LAPW method with a simplified version of the 
generalized gradient approximation [|18|] for the exchange-correlation functional in the LSDA 
description. This approximation takes into account the charge and spin inhomogeneity in a 
material by including gradient corrections in the energy functional and gives more precise 



ground state energy ||T9[ but does not improve the quasiparticle spectra JJO . 

The Figure 2 shows the calculated electron energy dispersions along the high-symmetry 
lines of the Brillouin zone (BZ). The most interesting feature of these calculations is a rela- 
tively large dispersion along the k z direction parallel to the Cu-0 chains and dispersionless 
bands perpendicular to them, i.e. along the k x and k y directions. Along the T-Z direction of 
the BZ, the two bands, which cross the Fermi level, split slightly, indicating that the inter- 
action between the two neighboring chains in the ladder is small, but still exists. The total 
density of electronic states of SrCu02 is given in Fig. 3 in the energy range between -7 eV 
and 10 eV. The electron density of states is calculated with 180 k points in the irreducible 



part of the BZ and it is in an agreement with previous publications |L0, yj. Because the 
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main properties depend on the electrons at the Fermi level, we paid a special attention to the 
energy dispersions near Ep. In the inset of Fig. 3 we show the partial density of electronic 
states in the energy range from -1.5 eV to 1 eV. Figure 4 gives the contribution from the 
different Cu 3d and 02p orbitals to the partial density of states in the energy range close to 
the Fermi level. From the results given in Figs. 2-4 we concluded that: 

(i) The electronic structure of SrCu02 is a quasi-one- dimensional one. 

(ii) The main contribution to the density of states in SrCuC>2 near the Fermi level comes 
from the Cu 3d y 2_ z 2 orbitals. 

(in) There is a small hybridization between the Cu 3d y 2_ z 2 and the O 2p z ,2p y orbitals. 

These results are in an agreement with the X-ray absorption spectroscopy data |2l| , which 
show that the holes in SrCu02 have predominantly the 3d y 2_ z 2 character (which is analogous 
to the Cu 3d x 2_ y 2 orbitals in notation for the 2D high-T c cuprates). The hole occupancy of 



3d 3x 2_ r 2 is less than 5% [|21|1 . 

However, the metallic state obtained in the framework of the LSDA calculation (in Figs. 2- 
4 the Fermi level Ep=0 is inside the occupied band) does not allow us to compare these 
results with our experiments. 

The applied LAPW method with the gradient LSDA corrections (or other LDA versions 
to the density functional) employs the orbitally independent exchange-correlation potential, 
which can not recognize different <i-orbitals in the copper ion with an open d-atomic shell, 
Cu 1+n (3d 10 ~ n ), of the SrCu02 material. As a result it gives a satisfactory band structure, 



but can not overcome " an energy gap problem" J22[ . In general, the local approximations to 



the density functional theory (for a review see and Refs. |23|, [2J, |25[ for an inclusion 
of the spin-dependent exchange-correlation potentials) are successful to describe the ground 
state energies and the quantities, weakly depending on a charge or spin density |2|, |7|, but 
the excited states are beyond the scale of the density functional theory. 

To alleviate an energy gap problem of L(S)DA to density functional theory the LDA+U 
method has been proposed (for details see Ref. |28j and its generalization, LSDA+U [f29fl). 
These methods have extended LDA-like functionals with an addition to them the on-site 
Hubbard terms and they are based on the unrestricted Hartree-Fock wave-functions. For 
the Mott-Hubbard insulating compounds with transition or rare-earth metals the orbital- 
dependent potential U leads to splitting of the partially filled 3d- or 4/-bands forming 
the upper and the lower Hubbard-like bands. The LDA+U functional gives more correct 
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antiferromagnetic properties than the LSDA one. Also LDA+U's are relevant for study of 
the temperature induced phase transition from an non-magnetic insulator to a metal with 
a local magnetic moments and for study of the charge orders, where LSDA fails. However, 
the used mean-field correction for a single orbital, Un^n^ — ► Un^n^in^ is an expectation 
value for the number of electrons with a spin projection j occupying a fixed partially filled 
orbital on a particular lattice site), to an exchange-correlation potential reproduces only 
qualitatively the correct physics for the Mott-Hubbard insulators. Inclusion of the inter- 
orbital intra-atomic interactions (Coulomb repulsion and Hund exchange) and often the 
antiferromagnetic ordering f^Oj (regardless the absence of an antiferromagnetism in material 
under consideration; e.g, SrCu02 is in paramagnetic state above 1.4K) does not change 
significantly the qualitative character of the LDA+U results. The obtained gap values differ 



by about leV from the experimental ones [28, 29 1. Also in the framework of LDA+U methods 



the band structure is unsatisfactory and these methods produce unphysical insulating states 
for transition metals since the LDA+U's split even the partially-filled bands due to the self- 
interaction corrections both for localized and extended electron states. It is also known that 
the electron correlations produce narrowing of the dispersion branches and consequently an 
increase of the effective masses [ftl , |3l^ , |33|| which is outside the LDA+U domain. 



As applied to SrCu02, the band gap values obtained in the LDA+U calculations [|TT|, [T^j 
differ significantly from our experimental findings (see Sec. IV). Besides that, a shift of the 
energy bands due to a finite U can not guarantee the correct band extrema positions to 



identify the direct or indirect character of the interband electronic transitions. In Ref. JTT 
the reported band gap 1.63 eV (the most close to our measured value) corresponds to an 
unrealistically large Anderson-Hubbard parameter £7=12.3 eV. 

Having all above mentioned in mind, we developed a model which introduces the electron 
correlations already in the zeroth order of the applied perturbation theory with respect to the 
hopping energies in a realistic lattice structure of SrCu02 (Sec. III. B). 

B. The tight binding method for the correlated electrons 



According to our LSDA data (Figs. 2-4) and also Ref. the main contribution to the 
electron density of states close to the Fermi level comes from the d y 2_ z 2 orbitals, slightly 
hybridized with the O 2p orbitals. In addition, our absorption measurements (see, Sec. IV) 
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show that the lowest energy gap is practically temperature independent (a charge transfer 
gap depends on temperature in accordance with the variation of the Cu-0 distance with 
temperature). All these important observations have lead us to the conclusion that the 
band gap in SrCuC>2 is a consequence of strong electron correlations, favoring the Mott- 
Hubbard like insulating state. 

Although the density functional theory reduces the Schrodinger equation for electrons and 
atomic nuclei in solids to non-linear single-particle equations (equations with self-consistent 
potential, which itself depends on the solutions of the equations), the explicit dependence of 
the ground state energy on the ground state charge and spin density is unknown and such an 
explicit functional may even not exist. As the starting many body linear Schrodinger equa- 
tion, the problem of the non-linear single-particle equations with the exchange-correlation 
potentials is a formidable problem which can not be solved without additional approxima- 
tions. The exchange and correlations are often approximated by the L(S)DA according to 
which the charge (spin) density in the exchange-correlation potential of an electron gas or 
a jellium model is replaced by the local density in a real material. The strategy of the 
LDA approach is close to the Landau fermi-liquid theory, where the energy is a functional 
with respect to a single-particle partition function. But contrary to the Landau theory, the 
density functional theory can calculate only the ground state energy, distributions of charge 
and spin densities and the quantities connected with them. The main problem in our case 
is that LDA and its versions (e.g., the one used by us) do not give a clear knowledge about 
the Anderson-Hubbard parameter (an essential on-site repulsion of the correlated electrons) 
in narrow-energy band materials such as SrCu02- 

The conventional tight-binding method (U = 0) for an electronic structure in solids is 
self-contradictory. According to it the eigen-functions of the Schrodinger equation are ap- 
proximated using the electron wave-functions of the isolated atoms, i.e. a conventional tight- 
binding method is the more suitable the greater interatomic distances in the crystal are. But 
in this case the prevailing terms in the Hamiltonian, the strong electron-electron interac- 
tions U, become even larger; they can not be reduced to any mean field version and then 
the problem is fully outside the domain of the standard Slater-Koster scheme. The present 
tight-binding method for correlated electrons is based on the hypothesis that the narrow- 
energy band material properties, e.g., SrCu0 2 in our case, are governed by the intraatomic 
electron correlations U , exceeding considerably the transfer energies. 
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A reasonable simulation of the many-body effects leads to tremendous problems in terms 
of the conventional Fermi- or Bose-operators permutation relations which do not encompass 
all possibilities of the second quantization formalism. The permutation relations for the 
Okubo-Hubbard X-operators are linear with respect to theirselves operators. The necessity 
to introduce the operators with more complicated permutation relations than the fermionic 
and bosonic ones for study of the electron correlations was indicated by Bogoljubov already in 
1949 R34 ]. We consider the tunneling part of any correlated Hamiltonian as a perturbation 
with respect to the strong electron correlations included in the unperturbed part of the 
Hamiltonian. The Hamiltonians with correlated electrons are rewritten in terms of the basis 
and only basis vectors of the corresponding superalgebra. To treatment of such Hamiltonians 
in the X-operator framework is based on a rigorous successive method which we followed 



in our previous publications |3^, [36|, [37JI an d it will be also applied here. The systematic 
perturbation theory is based on the generalized Wick's theorem as an iteration procedure 
reducing the time-ordered product of n of X-operators to product of n — 1 of thereof. The 
first order self-energy is the tunneling matrix itself from the perturbative Hamiltonian. 

Here in the framework of the su (2, 2) superalgebra approach for the SrCuC>2 system 
we will neglect the scattering of the correlated copper electrons by the spin and charge 
fluctuations, aiming at comparing the electron spectra with the conventional tight-binding 
calculations, which are done in the first order of the transfer energy. Contrary to the 
standard Hubbard model for the s-electrons with a single site per a unit cell, the starting 
Hamiltonian includes a realistic unit cell (a few sites) and the correlated electrons with 
the non-zero angular momenta. In the framework of the presented tight-binding method for 
the correlated electrons the non-spherical wave-functions provide anisotropy of the hopping 
integrals in the lattice. In the considered order of the perturbation theory we will concentrate 
on the influence of the band structure effects which are of significance for multicomponent 
systems such as SrCu02- 

The LSD A calculation of the electronic structure of SrCu02, Figs. 2-4, revealed that 
electron energy dispersions are governed mainly by electrons in the zigzag ladder. The 
interladder coupling is negligible because of the large (6.75 A) Cu atom distance between 
neighboring ladders. A unit cell for the correlated electron structure calculation includes 
two Cu-ions only, as indicated by a dashed rectangular in Fig. 1(a). We assume that a 
ladder unit Cu 1+n 02~ has the total charge —2, i.e. there is one hole, n — 1, per a copper ion 
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in the ladder. In our minimal model, we neglected hybridization with the p y ,p z — orbitals of 
the intermediate oxygen atoms, assuming the d y 2_ z 2 character of holes. Such approximation 



is common for the ID cuprates |38| because the Anderson-Hubbard repulsion usually opens 
a gap between the 3d bands. Inside the ladder, the arrangement of the Cu-atoms is such 
that the directions from the Cu-ion in one leg to the two nearest Cu-ions of the neighboring 
leg are almost at a right angle (~ 90°), see Fig.l. In our consideration of the electronic 
structure we assume that these directions form an ideal right angle and the electron energy 
dispersions are governed mainly by the correlated electrons in the single zigzag ladder. 

Aiming at the determination of the Anderson-Hubbard parameter U from the optically 
measured correlation gap values, the further theoretical study is based on the realistic Hamil- 
tonian with one copper rung/ dimer, a — b, per a unit cell: 
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where a, b denote chains (legs in the ladder directed along the c-axis), t is an amplitude 
of the carrier hopping along the legs, t yz is a diagonal hopping amplitude between the 
legs, U is the Anderson-Hubbard repulsion of the £^2_ z 2 -electrons on copper site and /i is 
the chemical potential. In LDA band structure calculations the electron-electron interac- 
tions are approximated by an effective single-electron problem, whereas in our approach 
the Anderson-Hubbard parameter U is included explicitely in the non-perturbative Green's 
function. Applying the A— operator machinery P7[ , the correlated energy bands are gov- 
erned by zeros of the inverse Green's function for the A-operators in the first order with 
respect to the tunneling matrix: 
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where r = — 2tcosp z , and v = — t yz (1 + e~* Pz ). Here the correlation factors, /o+(-2) , in 
the diagonal Green's functions are determined by the fermion occupation n per site. For a 
considered nonmagnetic phase of SrCu02 they are /o+ = 1 — n/2, f-2 = n/2 and all equal 
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to 1/2 (n = 1). After an analytical continuation, iu n — > £ + i8, from the secular equation 



D- 1 (u 



v 



one can find the four branches of the correlation energy dispersions in an 



explicit form: 



V ,2 



g(P) = 4 ,2 -yW+(x) (4) 

where e 1 ; 2 = —tcosp z ±t yz cos^. For the derivation of these energy dispersions from 
the four- by-four fold secular equation (see Eq.(2)) it was useful to apply the proved theorem 
about the decomposition of the determinant with respect to diagonal elements (see Appendix 
AinRef. ||). 

The lower correlated subbands, an d Ci> are completely occupied by the two holes from 
a unit cell of the zigzag ladder. It is essential that the chemical potential /i should be found 
self-consistently from the equation for a particle density 

n = T J2 e^D^p&Un), 

p,u},a,/3 

where a and f3 label the rows and columns of the Green's function elements (c.f. Eq.(2)). For 
n = 1 (the case of SrCuO^) the calculated chemical potential is positioned in the middle of 
the energy gap (/i=0). For other densities, n ^ 1, the used procedure generates a correlated 
metal. The upper correlated subbands, £ B , are empty for SrCuO^. The nearest unoccupied 
energy band is and the correlation gap (indirect transition, see Fig. 5) in an electronic 
structure can be estimated as 

A = min^ (p) - max^j (p) = 



Wu 2 +(|) + 



For the non-correlated energies, e 1,2 = 2e 1 ^ 2 (c.f., Eqs.(3,4)), the electron density of states 
per atom, po(e)= J2[${ £ ~ e^)+5(e — £^)}, is defined analytically as follows 

Vz 
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P ° iB)= r Jtl-U(e-2t)hjT=7f <6) 
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where 



si,2 = ^ (~fy* ± yjtf z -At (e-2t)) = -ar 4 ,3 



(7) 



The electron-electron repulsion, [/, splits the density of the non-correlated electron states 
Po(e). After the transformation to the "correlated" variables, Eqs. (3) and (4), one can get 
the correlated electron density of states as 
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Figure 5(a) shows the correlated electron energy dispersions (Eqs. (3) and (4)) along the Z 
and the Y symmetry directions of the BZ. The correlated electron density of states is given 
in Fig. 5(b). The overlap of the energy ranges for the energy dispersions, Eqs. (3) and (4), 
leads to the appearance of singularities at 
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and a takes sign +1 and —1 for the correlated subbands B (Eq. (3)) and A (Eq. (4)), respec- 
tively. The van Hove square root divergencies inside the correlated bands are manifestations 
of the quasi-one-dimensional electronic structure. 

The comparison of the calculated electronic structure with the measured optical transi- 
tions allows us to estimate the Anderson-Hubbard parameter and the hopping energies of 
correlated electrons in the zigzag ladder compound SrCu0 2 . 
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IV. EXPERIMENTAL RESULTS 



Figure 6 shows the pseudo-dielectric function e(uj) = e±(u) + ie 2 (uj) for a light polarized 
parallel to the a- axis (Fig. 6(a)), and the c-axis (Fig. 6(b)). As it can be seen from Fig. 
6(a), the maximum of e 2 (uj) is at about 1.77 eV. This value corresponds to the energy gap 
along the a-axis direction, which is perpendicular to the CuO chains. Besides the most 
intensive peak at 1.77 eV, we found the next peak at about 2 eV. To determine the precise 
energies of these electronic transitions we fitted both the real and the imaginary parts of 
the second derivative spectra d 2 e/dE 2 simultaneously by a least squares routine in terms of 
standard line shapes: 



where A, E and T are the amplitude, energy and half linewidth of the electron transitions, 
respectively. The is a phase factor and the exponent m has the value —1/2 for the one- 
dimensional case. The results of the calculations are given in the inset of Fig. 6(a). At 
higher energies we observed two additional peaks at about 4.1 eV and 5.1 eV. The origin of 
all these transitions will be discussed later. For E| |c polarization the first clearly pronounced 
peak appears at about 3.2 eV, and the next maxima is at about 5.2 eV (see Fig. 6(b)). For 
this polarization the low intensity of the dielectric function does not allow us to apply the 
same fitting procedure as in the E||a case. 



An unpolarized reflectivity spectrum is shown in Fig. 7(a) in the energy range from 
1 to 5.5 eV. The peak positions at 1.42 eV, 1.73 eV, 2.04 eV, 2.5 eV, 2.9 eV, 4 eV and 
5.2 eV are determined as the maxima of an optical conductivity function obtained from 
the Kramers-Kronig analysis of the reflectivity spectrum. In order to distinguish different 
polarization contributions in the reflectivity spectrum shown in Fig. 7(a), we calculated the 
polarized reflectivity spectra using the pseudo-dielectric function data from Fig. 6. The 
results of the calculations are given in Fig. 7(b). As it can be seen from Fig. 7(b), the main 
contribution to the unpolarized spectra comes from the E||a polarization. As a proof that 
the lowest intensity peak in unpolarized reflectivity spectra belongs to the E||c polarization, 
we measured polarized photoreflectance spectra in the energy range below 1.5 eV, which 
are shown in the inset of Fig. 7(a). These spectra clearly demonstrate that the reflectivity 
maximum at the lowest energy in Fig. 7(a) originates from the E||c polarization. The peak 
position of this transition is found at 1.43 eV using a photoreflectance fitting line shape 



e (u) = C - Ae^iuj - E + iT) m , 
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procedure based on the third derivative of the primary spectra |40| . Thus, we conclude that 



the lowest energy gap in SrCu02 is at about 1.42 eV for the polarization along the Cu-0 
chains. 

Figure 8 represents the absorption spectra of SrCuC>2 measured at room temperature 
and 5K for the incident light polarized along the c- and the a-axis. These spectra were 
calculated using the relation a(u) = (l/d)ln[(l — R(u)) 2 /T(u>)}, where R and T represent 
the reflectivity and transmission coefficients, while d is the thickness of the sample. As it 
can be seen from the Fig. 8, there is a strong anisotropy in the positions of the absorption 
edges for these polarizations. 

V. DISCUSSION 

The experimental results given in Figs. 6-8 can be summarized as follows: 

(i) A strong anisotropy of optical properties is clearly seen. In the direction perpendicular 
to the Cu-0 chains (E| |a) the spectral weight is centered at about 1.8 eV. All other transitions 
at higher energies have smaller contributions to the dielectric function for this polarization. 
In the case of the E||c polarization (Fig. 6(b)) the spectra show very low intensity and the 
spectral weight for this polarization is shifted to higher energies at about 3.2 eV and 5.2 eV. 

(ii) The lowest energy gap appears at about 1.42 eV, well bellow the charge transfer gap 
(1.7 - 2.0 eV) of other CuO based materials. 

(Hi) The lowest energy gap for E| |c shows a very small energy change as the temperature 
is lowered, whereas the absorption edge of the corresponding gap for the E| |a shifts to higher 
energies. 

(iv) The higher energy electron transitions appear at about 4.1 eV and 5.2 eV for both 
polarizations. 

These features can be successfully explained in the framework of the electronic structure 
calculations, Sect. III.B. The parameters of our model for correlated electrons, t, t yz and U, 
were fitted to the experimental value of the indirect A gap (A =1.42 eV) and an exchange 
energy along the chain (leg) J = 0.18 eV p. We used also the ratio of exchange energies 
J'/J=0.16 (J' represents the exchange energy between legs), because this value is in the 



middle of the range J'/ J = 0.1 — 0.2 proposed for this kind of copper oxides . Using 
the relation between exchange and hopping energies J = At 2 /U and Eq. (5) for the lowest 
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energy gap, we estimated the other model parameters, £=0.30 eV, £^=0.12 eV, U = 2.0 eV. 
The obtained magnitude of the Anderson-Hubbard parameter of SrCuC>2 somewhat lower 
then the one (U = 2.1 eV) obtained previously by us in the electronic structure study of 
Sr 14 Cu 2 40 4 i P5| . The hopping energy t=0.3 eV along the legs is larger than in the case 
of Sr 14 Cu 2 40 4 i (£=0.26 eV) due to the difference between the exchange coupling constants: 
J = 180 meV in SrCuC>2 is larger than the J = 128 rneV value in Sri 4 Cu240 4 i. The fact 
of the coincidence of the Anderson-Hubbard parameters U for SrCuC>2 and Sri 4 Cu240 4 i is 
natural to expect because the local environment of Cu 2+ ions is similar in these oxides and 
the on-site Coulomb repulsion is governed mainly by the electron interactions on the same 
copper (i-orbital. 

The energy vs. wavenumber curves for correlated electrons (Fig. 5 (a)) show a large 
dispersion along the k z direction (parallel to the Cu-0 chains, along the c-axis) and no 
dispersion along the k x direction, in agreement with the LSDA calculations (Fig.2). Never- 
theless, the lowest energy gap at about 1.42 eV appears for the transition from the Z to the 
T point of the BZ (indirect transition, denoted as Ao in Fig. 5 (a)). The direct transition 
gaps A 2 , A 3 and A 4 at T (X) point of the BZ, are found to be 1.88 eV, 2.05 eV, and 2.19 
eV, respectively (see Fig. 5(b)). 

According to our calculations in Sec. III. B, the correlation gap A , Eq.(5), represents the 
lowest energy transitions between the split 3d y 2_ z 2 — 3d y 2_ z 2 states within each Cu atom. 
Thus, it is natural to expect a negligible temperature dependence of the correlation gap. In 
fact, experimental data presented in Fig. 8 show no temperature shift of the absorption edge 
for the E||c polarization. It means that SrCu02 belongs to the group of low-dimensional 
insulators with a correlation gap. Further support for this assumption can be found in 



the case of S^CuOs, also 1 D cuprate. Maiti e£ al. [f|2| assigned the gap, observed at 
about 1.5 eV in Sr 2 Cu03, as an insulating gap of the correlated nature. Nevertheless, our 
Raman spectroscopy study clearly demonstrated that the resonance behavior in SrCu0 2 
differs strongly from that in the 2D cuprates. Namely, the resonant Raman scattering shows 
that the infrared active modes and their overtones resonate much more strongly for the 
laser energies near the charge transfer-gap in the 2D insulating cuprates ||. In ID SrCu0 2 
the infrared modes resonate at energies noticeably higher than the correlation gap A 0. 
Finally, we have found that a similar spectral weight distribution as the one given in Fig. 
6(b) is also reported for Li 2 Cu02 for the light polarized along the chains [fi5|] . This supports 
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once again our findings that the electronic structure of the ID cuprates differs significantly 
from that of the 2D cuprates. 

In the inset of the Fig. 5(a) we show the A gap dependence on the U parameter, Eq.(5), 
for the two different values of the interchain hopping t yz at £=0.30 eV. The limiting case 
t yz =0 represents one leg of the ladder or a single-chain structure which exists in Sr 2 Cu0 3 . 
As it is seen in the inset of Fig. 5(a), the interchain hopping reduces a gap. It means that 
the corresponding gaps in the single-chain compound Sr 2 Cu03 should be somewhat larger 
than in SrCu0 2 . This agrees well with our experimental findings. Namely, the absorption 



edge in Sr 2 Cu03 is at 1.5 eV J42|, about 5% higher than in SrCu0 2 , see Fig. 8. 

The electron structure calculations for correlated electrons (Fig. 5) make it possible to 
compare to experimental data not only the lowest energy gap value. The peaks at 1.77 
eV and 2 eV of the E||a £ 2 (u;) spectra correspond to the L 2 „-Li + and L 2 _-L 2+ transitions 
at T and X-points (A 2 and A 3 gaps), respectively. In the E||c polarized £ 2 (u;) spectra we 
observed three very low intensity peaks at about 2.2 eV, 2.6 eV and 3 eV. These transitions 
can be assigned as L]_— L i+ (A 4 ), L 2 _— L 3+ and Li_— L 3+ , respectively. Because of the very 
low intensity of these modes, it was hardly possible to extract their exact positions from the 
noise level. The energies of all possible electron transitions between occupied and empty 
states are shown in the inset of Fig. 5(b). 

Now we discuss the two higher energy features at about 4.1 eV and 5.2 eV. They have 
been already observed in a lot of cuprates, as discussed in Ref. |44|| . Alonso et al. HI 



shown that in Nd 2 _ a; Ce :1 .Cu04 the transition at about 4 eV originates from the Cu 3d - 
Cu Ap intraionic transition and the highest energy transition about 5.2 eV occurs due to 
transition between the O 2p and the Cu 4s states. The existing band structure calculations 
of SrCu0 2 are not detailed enough to corroborate this assignment. 

Let us now consider Fig. 8, where the absorption coefficient spectra are presented. There 
is a large difference between the absorption edge positions for the E| |c and the E| |a polariza- 
tions. Besides that, the slope of the absorption coefficient is higher for the E| |c then the E| |a 
polarization, as a consequence of a quasi ID character of this transition. A general approach 
to analyze the fundamental absorption edge is based on the use of the power-law behavior 
of a{uj) in the vicinity of the band gap, E gap . The absorption coefficient can be described as 
a(uj) = A(uj — E gap ) k , where A is a slowly varying function which is regarded as a constant 
over the narrow range under study, and a number k depends on the nature of an electron 
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transition from the occupied to the empty band. The k value is 1/2 for a direct transition 



and 2 for an indirect one. Thus, the ya(u) ~ (u — E gap ) dependence can be used to obtain 
the indirect gap value by extrapolating the linear portion of this curve to intersect with the 
uj axis. The direct gap position can be extracted from a maximum of the first derivative of 



a(tu) (a'(uj) ~ 1/y^ — E gap ). This procedure is valid for a three dimensional case and for 
the parabolic electron energy dispersions. 

The absorption spectra for strongly correlated electron systems were analyzed using the 



density of states vs. energy dependence of ID [45 or 3D [[|(| semiconductors with parabolic 
energy dispersions. Actually, the dispersion branches in strongly correlated systems are 
parabolic in the vicinity of the high symmetry points of the BZ: cos p = 1 + p 2 /2 around T 
(~ 0) and cos p = l—p 2 /2 around Z (~ 7r). Furthermore, from the extrapolation of the linear 
part of the a(u) to the energy axis, we obtain the value of 1.52 eV (1.65 eV) for E||c (E||a) 
polarization. This energy is higher (lower), than that obtained by ellipsometric or reflectivity 
measurements. The complete agreement with ellipsometric and reflectivity data is achieved 
when we consider the E\ \ c absorption edge as an indirect transition and the E\\a as a direct 
one, which is illustrated in the left and the right insets of Fig. 8. Thus, we have concluded 
that the correlation gap in SrCu02 represents the indirect transition of carriers between the 
occupied and the empty correlated subbands. The indirect transition requires a change in 
both energy and momentum of carriers, as our band structure calculations predict, see Fig. 
5(a). It is well known that in semiconductors an electron momentum is conserved via an 



interaction with phonons p6| . For such conclusion here further experiments are necessary to 
clarify this point. They can be done by the transmission measurements through the samples 
of different thickness and/or comparison of photoluminescence with photoreflectance spectra. 

Quite recently, based on the Bethe ansatz solution for the Hubbard chain it was shown 
|47| , that in an insulating state the optical conductivity can be described as o-(uj) ~ C(u — 



Emh) 1 ^ 2 , where Emh is the (Mott-Hubbard) gap. Since a(cu) ~ cr(u>), this dependence has 
the same power-law as the one we used for the direct electron transition. 

For the E||a polarization, the absorption edge shifts to higher energies by about 0.06 
eV, when the temperature decreases from a room temperature to 5 K. For this electronic 
transition we cannot neglect hybridization of the Cu 3d with the O 2p bands. Thus, the 
absorption edge shifts to higher energies according to the change of the Cu-0 distance 
with temperature. Besides that, this energy gap represents the transition between the 
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dispersionless branches with very high effective electron masses and, consequently, with 
a strong influence of ligands. Because of that we believe that for the E||a polarization the 
absorption edge shift appears mainly due to dilatation of lattice by lowering the temperature. 

In conclusion, we have investigated the electronic structure of the zigzag chain compound 
SrCuC>2 combining polarized optical absorption, reflection, photoreflectance and pseudo- 
dielectric function measurements with theoretical estimations. At 300 K these measurements 
yield the energy gaps 1.42 eV and 1.77 eV along and perpendicular to the Cu-0 chains, re- 
spectively. The electronic structure of this oxide is calculated using the LSDA with gradient 
corrections and the tight-binding method for the correlated electrons. The gap value of 
1.42 eV (1.86 eV) is found for the electron hopping energies between copper sites along 
legs, t = 0.30 eV, and between them, t yz — 0.12 eV, with the Anderson-Hubbard parameter 
U = 2.0 eV. The obtained experimental results and electronic structure calculations have 
shown that the SrCu0 2 zigzag chain compound belongs to the low- dimensional insulators 
with the band gap of the correlated nature. 
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FIG. 1: Crystal structure of SrCu02 in the (a) (100) and (b) (001) plane. Dashed line rectangle 
represents a unit cell of the zigzag ladder. 

FIG. 2: Electronic structure of SrCu02 along the several high symmetry directions in BZ calculated 
using the LSDA with gradient corrections method. The cartesian coordinates of high-symmetry 
points are as follows: T (0,0,0); Y(2vr/a,0,0), (0,2tt/6,0); Z(0,0,tt/c); T(27r/a,0,vr/c), (0,2vr/6,vr/c). 

FIG. 3: Total density of states calculated within the LSDA with gradient corrections method. 
Inset: Partial density of states close to Fermi level of different atoms of SrCu02- 

FIG. 4: The contribution of different orbitals to the calculated partial density of states (DOS) of 
SrCu0 2 . (a) Cu-d DOS, (b) Ol-p DOS, and (c) 02-p DOS. 

FIG. 5: (a) The tight-binding dispersions for correlated electrons in SrCu02 with parameters t=0.3 
eV, t yz =0.12 eV, U=2.0 eV. The momenta are given in units | p y V% |=| p z \=n of the Brillouin 
zone boundaries, the Fermi energy Ep = is inside of the correlation gap. Inset: energy gap vs. U 
dependence for t yz = and 0.12 eV; (b) The correlated electron density of states as a function of 
energy. The electronic structure parameters are the same as for Fig.5(a). Inset: the energy values 
of all possible transitions between occupied and empty states. 

FIG. 6: Room temperature real (ei) and imaginary (e 2 ) part of the pseudodielectric function of 
SrCu02- The spectra of the (010) surface taken with (a) the a-axis (E||a) and (b) the c-axis 
(E||c), parallel to the plane of incidence. Inset: second derivative of dielectric functions for E||a 
polarization in the 1.5 - 2.5 eV spectral range. 

FIG. 7: (a) Room temperature unpolarized reflectivity spectra of the SrCu02 single crystal, (b) 
Reflectivity spectra calculated using experimental data of the pseudo-dielectric function from Fig. 
6. Inset: polarized photoreflectance spectra at room temperature. 

FIG. 8: Absorption coefficient spectra of SrCu02 at room and liquid helium temperatures. Inset, 
left panel: ^a(oj) vs. uj dependence. Inset, right panel: first derivative of a(uo) dependence for 
E||a polarization. 
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